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1  Introduction  and  Summary 


During  the  first  year  of  this  program,  we  focussed  on  development  of  an¬ 
alytical  methods  for  representation  of  multiple  scattering  effects  in  het¬ 
erogeneous  media.  We  have  examined  several  methods  for  representation 
and  computation  of  effective  parameter  models  for  the  media  in  terms  of 
statistical  characterizations  of  the  media.  These  methods  lead  to  represen¬ 
tation  for  the  effective  parameters,  e.g.,  the  effective  dielectric  constant  or 
conductivity  in  a  heterogeneous  medium.  Since  these  representations  may 
be  difficult  to  evaluate  in  specific  cases,  we  have  also  examined  current  re¬ 
search  on  the  derivation  of  bounds  for  the  effective  parameters  as  a  means 
of  quickly  providing  good  approximations  to  the  behavior  of  the  medium. 

Our  work  on  the  T-matrix  methods  is  summarized  in  section  1.  We  plan 
to  continue  this  work  to  develop  numerical  algorithms  for  the  evaluation  of 
scattering  cross  sections  and  absorption/dissipation  mechanisms  in  typical 
foliage  configurations.  It  is  important  to  emphasize  that  the  T-matrix 
methods,  though  formal  in  structure,  provide  simple  models  for  multiple 
scattering  phenomena  which  are  valid  for  a  wide  range  of  frequencies  (of 
the  incident  radiation). 

The  work  on  the  use  of  “multiple-scale"  asymptotic  analysis  (homoge¬ 
nization]  for  the  representation  of  scattering  processes  is  given  in  section 
2.  This  is  primarily  a  “low  frequency”  theory,  valid  when  the  wave  length 
of  the  incident  radiation  exceeds  the  characteristic  dimensions  of  inter- 


a 
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particle  aeparation  in  the  scattering  medium.  The  homogenization  method 
does,  however,  provide  a  systematic  procedure  for  the  construction  of  high 
order  approximations  to  the  “effective  parameter”  representations  of  the 


scattering  process.  These  approximations  may  be  evaluated  by  common 


numerical  procedures  when  the  scattering  medium  has  a  periodic  struc¬ 
ture.  We  describe  this  case  fir^t  in  section  3,  before  turning  to  the  more 
realistic,  but  less  tractable,  case  when  the  physical  model  parameters  (per¬ 
mittivity,  conductivity,  and  permeability)  are  random  functions  of  position 
in  the  medium.  We  have  not  completed  our  analysis  of  the  probabilistic 
case  at  this  writing,  and  we  report  only  on  a  basic  mathematical  result  for 
a  key  quantity  arising  in  the  homogenization  method. 


In  section  4  we  examine  the  derivation  and  use  of  “bounds”  on  the 


effective  parameters  in  a  scattering  problem.  The  derivation  of  bounds 
is  an  old  subject  in  scattering  theory,  dating  from  early  in  this  century. 
Nevertheless,  there  has  been  some  very  interesting  recent  work  on  various 
methods  for  deriving  bounds.  We  discuss  this  work  briefly  in  section  4,  as  it 
applies  to  the  problem  of  scattering  from  foliage.  As  with  homogenization, 
the  class  of  bounds  disciissed  applies  mostly  to  the  low  frequency  case. 


2  T-matrix  Approximations  to  Multiple  Scat¬ 


tering 

In  this  section  we  consider  a  classical  approach  based  on  the  “T-matrix 
formalism”  for  the  problem  of  modeling  scattering  in  a  region  containing 
two  or  more  types  of  dielectric  scatterers  arranged  in  a  random  geometry 
in  a  surrounding  medi\im.  Models  for  scattering  from  a  medium  containing 
a  two  (or  more)  classes  of  (randomly  oriented)  scatterers  in  an  enveloping 
medium  are  prescribed  in  terms  of  the  definition  and  derivation  of  approx¬ 
imations  for  the  “effective  dielectric  constant"  resulting  in  a  family  of  ap¬ 
proximations  for  an  “effective  scattering  representation”  for  heterogenous 
media. 

2.1  Basic  Framework 

Consider  a  region  0  of  space  occupied  by  a  homogeneous  (background) 
dielectric  medium  with  dielectric  constant  «o  containing  elements  from  two 
classes  of  scatterers  with  dielectric  constants  ci,£j.  Let  t;„,x  =  l,2,j  = 
1,2, .. .  ,iVi,  be  the  subsets  of  0  occupied  by  scattering  elements  j  of  class 
i.  The  dielectric  properties  of  the  composite  material  in  0  are  described 

by  ^ 

c(*)  =  Co  +  51 (c,  -  £o)  (*)  (i) 

i=i  j=i 
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aTjC, 


Xw,(3:)  = 


1  I  fc  Vij 

0  X  ^  Vij 


Let  6i  be  a  (dimensionless)  parameter  describing  a  characteristic  dimension 
of  elementary  scatterers  in  class  i  (e.g.,  the  radius  of  spherical  scatterers); 
and  let  p,  be  the  total  volume  fraction  of  0  occupied  by  scatterers  of  class 
t  =  1, 2.  Suppose  a  constant  (in  amplitude  -  time  harmonic)  field  is  incident 
on  the  region  0. 

We  wish  to  characterize  the  scattering  properties  of  the  composite  ma¬ 
terial  in  the  limit  JVj  -♦  oo,  6i  — ►  0,*  =  1,2,  with  p\,P2  constant.  We  are 
particularly  interested  in  evaluating  the  relative  roles  of  the  two  classes  of 
scatterers  in  the  scattering  process,  e.g.,  in  terms  of  their  relative  densities 
Pi  in  the  cases  (a)  pi,P2  —*  0  with  pi/pj  ^  0(1)  and  (b)  pi,P2  0  with 

PilPj-^0. 

We  shall  sketch  the  development  of  a  family  of  effective  scattering  ap¬ 
proximations  for  the  material  using  the  multiple  scattering  framework  de¬ 
veloped  originally  by  M.  Lax  (1951,1952). 


2.2  Transition  Operator  Representations 

We  use  the  transition  operator  formalism  of  classical  scattering  theory  (Lax 
1951,1952,1973)  as  used  by  (Lang  1981)  and  (Kohler  and  Papanicolaou 
1981)  as  a  framework  for  the  derivation  of  several  formulas  for  an  effective 


dielectric  constant  for  the  composite  medium.  Maxwell’s  equations  for  this 
situation  may  be  written  in  the  form 


V  .  (eoE)  +  -  £o)  E  V  •  [x.,,E]  =  0  (3) 

.=i  ,=i 

V  X  ^  =  0 


with  the  boundary  condition 


<  E  >—*  E  as  Ni,N2—*oo  (4) 

where  <  •  >  is  expectation  and  E  is  the  constant  external  field.  The 
constituent  relationship 

D(i)  =  <(i)£(i)  (5) 

and  the  condition 

<  £>(•)  >=  e'E  (6) 

defines  the  effective  dielectric  constant  for  the  composite  material. 

Rewriting  (3)  in  abstract  form 


(Lo  +  M)E  =  0,  V  X  £;  =  0  (7) 

where 

Lo  =  V  •  (co‘) 

AT. 

Af  =  Afi  +  Mj,  Mi  =  Vii  (8) 

j=i 
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V«  =  (<i-<o)-Vlx.,,l 

we  have 

E  +  L^^ME  =  F  (9) 

as  the  integral  form  of  (3).  F  Is  chosen  so  that 

LqF  =  0,  <.  E  >=  E. 

The  transition  operator  (“matrix”)  T  is  defined  as  follows: 

=  (Lo  +  M  -  M)-^M 

=  [(Lo  +  Af)(/  -  (Xo  + 

=  r(/-r)-^ 

with  T  =  (Xo  +  So  (9)  becomes 

[/  +  r(/  -r)-‘]x  =  F 

or 

E  =  [I-T)F 

The  condition  (10)  gives 

F  =  (/-<r>)-‘x 

E  =  (I-T)(I-  <  T  >)-^E 
Using  (5)  (6)  (15),  we  have 


(•E  =(<€>-  <eT  >)(/-  <  T  >)-^E 


as  the  definition  of  «*  in  terms  of  the  transition  operator.  It  is  apparent 
that  knowledge  of  T  completely  characterizes  the  scattering  properties  of 
the  composite  medium.  The  operator  T  is  related  to  the  scattering  am¬ 
plitudes  of  the  individual  particles,  and  to  their  interactions.  For  simple 
(single)  dipole  scatterers,  T  can  be  determined  from  the  polarizability  of 
the  element;  and  so,  T  is  defined  in  terms  of  quantities  of  physical  in¬ 
terest.  Other  properties  of  the  transition  operator  we  described  in  (Laix 
1951,1952,1973),  (Lang  1981),  and  (Kohler  and  Papanicolaou  1981).  Our 
analysis  of  the  representation  (15)  (16)  is  based  on  (Kohler  and  Papanico¬ 
laou  1981). 

To  bring  out  the  roles  of  the  two  classes  of  scatterers,  we  reconsider  (9) 


E  +  4-  M2)E  =  F 


which  may  be  written  as 


E  =  F  -  TiEi  -  T2E2 


=  r.(/  -  r.)-‘ 

Ti  =  (Lo  +  M.)-*M, 


E,  =  {I-T\)-^E 


since  Ei  —  TiEi  =  E,  y/e  have 


Ei  +  j2TtE,  =  F 


(The  sum  in  (20)  is  degenerate;  but  this  is  the  appropriate  form  when  0 
contains  K  >  2  classes  of  scatterers.)  Equations  (18)-(20)  express  the  field 
in  terms  of  the  fields  and  transition  operators  associated  with  each  class  of 
scatterers  in  the  composite  medium. 

Now 

T={Lo-\-Mi+M2)-^{Mx  +  M2)  (21) 

=  {I  +  T,2)-^Ti  +  {I  +  T2i)-^T2 
where  the  pair  transition  operators  are 

Tij  =  (£o  +  ^  y  =  1, 2.  (22) 

This  representation  gives  the  formula 

c*  - 1<  €>-<€(/  +  ri2)-^ri  >-<€(/  +  ri2)-'r2  >]  (23) 

■[/-  <  (f  +  T,:)-'T,  >  -  <  (/  +  Tu)-'T!  >]-' 

for  the  effective  dielectric  constant  in  which  the  roles  of  each  class  of  scat¬ 
terers  amd  their  interaction  is  explicit.  We  could  continue  this  process  to 
make  the  role  of  each  elementary  scatterer  of  class  i  explicit.  Note  that 
(23)  contains  no  approximations. 

2.3  Average  T-Matrix  Approximation  (ATA) 

Expanding 

T  =  Ti  +  T2-  T12T1  -  T21T2  +  •  •  •  (24) 
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neglecting  the  higher  order  terma  {T12T1,.. .),  and  substituting  in  (23),  we 
obtain  the  simple  apprximation 

«•-[<«>-  ^  <  tr,  >\[I  -  ^  <  T,  >]-»  (25) 

>=i  «=i 

called  the  average  T-matrix  approximation  (ATA)  (Lax  1951,1973)  (Kohler 
and  Papanicolaou  1981).  The  expressions  <  Tj  >  and  <  eTj  >  in  (25) 
diverge  in  typical  cases,  so  the  ATA  must  be  used  with  care. 

It  is  possible  to  evaluate  (25)  in  the  limit  as  the  sizes  of  the  elemen¬ 
tary  scatterers  approaches  zero.  Suppose  each  class  consists  of  spherical 
scatterers  having  radii  ^,,i  =  1,2  with 

Pi  =  (26) 

the  volume  fraction  occupied  by  each  class.  Here  c,  is  the  average  number 
of  sphere  centers  per  unit  volume,  this  case  leads  to  the  approximation 


•  J  *  1.2co+<» 

F  ■  p  ■  ■  m  m  —M 


1 


<Q  ^<1 
.2fO‘*-C3  j 


which  is  a  version  of  the  Clauaaius  -  Moaaotti  formula  (Kohler  and  Papani¬ 
colaou  1981,  p.  213).  In  general,  the  Claussius  -  Most»^  ■  ^i  formula  is  viewed 
as  a  good  approximation  for  low  to  moderate  volume  fractions  p,-,i  =  1,2. 

The  ATA  not  only  contains  divergent  terms  as  a  rule,  it  also  fails  to 
account  for  the  interactions  between  scatterers.  A  family  of  more  refined 
approximations  has  been  developed  to  include  interaction  effects  in  formulas 
for  effective  scattering  representations  which  are  nonetheless  computation¬ 
ally  feasible  (Lax  1951,1952,1973)  (Elliott  et  al  1974). 


.rr.r. 


2.4  Coherent  Potential  Approximation 


The  simplest  of  these  is  the  coherent  potential  approximation  (CPA)  which 
results  (in  effect)  from  neglecting  the  difference  between  the  field  exciting 
the  medium  and  the  average  field.  Let  he  parameters  to  be  selected, 

and  consider  Maxwell’s  equation  rewritten  as 

x;  V  •  lAtoEi  +  E  i:  V  •  i(£i  - 

•=1  «ali=l 

+  EI:-S-(1-A)V-£  (28) 

«=i  i=i 

V  X  i?  =  0 

We  may  write  this  abstractly  as 

=  0  (29) 

with 

=  V  •  (/?iCo')  +  V  •  [0tea’) 

=  (30) 

.=1  j=i 

Vf,=V-  ((££  -  £o)x«-|  +  ^(1  - 


The  same  argument  as  before  yields  the  representation 

~  1<  ‘  >  -EE  <  >lli  -  EE  <  I)’  >!■■ 

•=l j=l  <=1 j=l 

where 


lif  + 


(31) 


(32) 
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The  CPA  is  based  on  choosing  0  to  “optimize”  the  approximation  in  (31). 


Let 

w‘>(z)  =  (TfisUx) 


SO 

V  .  +  V  .  ((e,  -  Co)x«.,ti;’’  +  ^(1  -  0iW’] 

=  V  •  [(£.-  -  co)x»„‘/  +  ^{1  -  0i)9] 

V  X  w'^  =  0 


Let 

ef'^  =  [0i€o  +  1^(1  -  /?i)]eo  «  /?<eo  for  Ni  large 
and  =  w'x  +  vjj  where 

V*(£f’^<)  +  V.[(£,-£o)x«.,«;ri 

=  V  •  ((€i  -  €o)x«.,<7l 
V  .  +  V  •  [(fii  -  to)Xv.,«^i1 

V  X  =  0  =  V  X 


(33) 

(34) 

(35) 

(36) 

(37) 


Let  {y,^}  be  the  locations  of  the  centers  of  the  individual  scattering 
elements.  For  small,  we  may  replace  g  on  the  right  in  (36)  by  g{y^)  =  yi>. 
Using  this  and  the  approximation 

£f’^  <55  /3j€o  for  Ni  large  (38) 

11 


we  have 


W 


a 


gp  -  c3t7  \sii  ‘  -  yjj) 

3y3i£o  +  £,-  -  €o  ‘  ‘  |a:  - 


Jpi£o  +  £*•  —  £o 

a  result  similar  to  (Kohler  and  Papanicolaou  1981,  sect.  8).  In  addition, 
xising  (38)  in  (37) 

where 

•i  /  \  €o  “  0i)  r3t^  )  \_  N\  ^  c  fAt\ 


Combining  these,  we  have 


(l-^.)„  r{x-y)-g{y) 


li-yl^ 


h  3/?<€o  +  g,  -  gp  ’  \x  -  y^f|3 


~  T.  L  X«o  (x)ff.-,  +  E  li  ^3  (®) 

i=i  '*H»gp  -r  g.  gp  ,_i  i_i 

Arguing  as  in  (Kohler  and  Papanicolaou  1981,  sect.  8),  the  last  term  on  the 
right  in  (42)  is  0(JV,“‘)  on  the  average  as  iVj  — ♦  oo.  We  take  the  expectation 
and  choose  0  so  that 

<'tE^fit>^o  («) 

•=i  ,=l 


Thb  is  the  CPA.  This  requires 

(1  -  A 


<0  -  c. 


Z0i  ^•  3Aco  +  £.- -Co’* 
where  p<  =  with  Cj  the  average  number  of  scattering  centers  (of  class 

t)  per  unit  volxime. 

With  0i  chosen  from  (44),  we  have  <  >=  0,  and  so, 

the  effective  dielectric  constant  may  be  approximated  using  the  CPA  by 
evaluating  (31). 


3  Homogenization  and  Multiple  Scattering 

In  this  section  we  turn  our  attention  to  an  alternative  method  for  represen¬ 
tation  of  multiple  scattering  effects  in  heterogeneous  media.  The  method  is 
a  variation  of  the  “homogenization”  procedure  which  has  been  used  widely 
in  mathematical  physics  and  engineering  to  develop  effective  media  approx¬ 
imations.  This  method  has  promise  for  for  several  reasons: 

1.  It  leads  more  naturally  than  the  T-operator  formalism  to  a  com¬ 
pact  representation  for  an  effective  parameter  (complex  dielectric  con¬ 
stant)  representation  for  the  effective  medium  approximation. 

2.  The  resulting  representation  for  the  effective  dielectric  constant  in¬ 
cludes  the  interaction  of  microscopic  effects  (multiple  scattering)  ex¬ 
plicitly. 


3.  The  underlying  analysis  ap;  lies  (however,  with  significant  differences 
in  detail)  to  both  periodic  and  random  media.  In  the  former  dase  it  is 
possible  to  solve  the  equations  for  the  approximation;  in  the  latter  it  is 
necessary  to  develop  approximations  which  use  physically  measurable 
quantities  (second  order  statistics  and  correlation  functions). 

4.  One  can  prove  convergence  of  the  scaled  model  to  the  “homogenized” 
model  in  both  the  periodic  and  random  case. 

5.  The  representation  provides  a  (formal)  basis  for  the  systematic  con¬ 
struction  of  a  sequential  approximation  to  the  effective  dielectric  con¬ 
stant,  including  “higher  order”  expansions  for  the  effective  parameters 
in  terms  of  the  small  parameter,  at  least  in  the  periodic  case. 

We  wish  to  undertake  a  systematic  investigation  of  such  approximations 
in  the  context  of  scattering  from  foliage  covered  terrain,  and  the  multiscale 
(homogenization)  method  offers  a  more  general  setting  for  such  a  compar¬ 
ison  than  does  the  T-matrix  formalism,  at  least  in  the  frequency  range  to 
which  it  applies. 

3.1  A  General  Model 

To  illiistrate  the  ideas,  consider  the  following  general  model:  Let  0  C  be 
a  region  in  which  e(x),/i(z)  and  a(x)  are,  respectively,  the  dielectric  tensor, 
magnetic  permeability  tensor,  and  conductivity  tensor  (3x3  matrices)  of 
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the  material  in  region  0.  In  0  the  (vector-valued)  electric  and  magnetic 


9 

L 

I 

•2 

I 


fields  satisfy 


d  e{x)  0  E{x) 

*  0  ^(x)  H{x) 


0  Vx  E{t,x) 

— Vx  0  H{t,x) 


a{x)  0  E{tyx) 


[0  oj  [fr(t,x)  J 

E{Q,x)  =  Eo{x),HiO,x)  =  Ho{x)  (46) 

V.(ci;ol=0,V.[Mf?o]=0  (47) 

where  Vx  is  the  curl  operator  and  V*  =  div.  Note  that  (47)  implies 

V.[€£;(t,i))  =0,V.(MH(t,x)l  =0,Vt  >0 

It  is  necessary  to  assume  that  c(x),m{x)  and  a(x)  are  symmetric,  and  that 
c(x),/i(x)  are  positive  definite  matrices. 


Using 


0  _ -2- 

^  Szt  Si, 


VxE=  ^ 

9z, 


5^  El 

9zi  * 


-2-  0 

9z,  Szi  ^ 


we  see  that  (45)  may  be  written  in  the  abstract  form 


where 


EltyX)  _ 

v{t,x)=  ^  €»• 

H{t,x) 


iirjMPyBPEpBWB^WBnWBHHWW^*BvfBfflBtBWgrCTSrCT^^^^^?rTBWCT5BBB5BBRtTvivCcHff6cCcffff^y 

*W!«yWcMMWvSS«www8BWWff .  •  »w\wsvvvwWSisS  w>Kw?*S»!?SS»'Sv3?»!w 


A\x)  = 


t{x)  0 


L  0  J 

the  ,  i)  are  skew-symmetric,  6x6  matrices  with  O’s  and  I’s  as  elements, 


B(x)  =  ^  ^  6 

0  0 


The  symmetry  assumption  means  that  (49)  is  a  symmetric  -  hyperbolic 
system  of  a  type  common  in  mathematical  physics.  The  (asymptotic) 
ansdysis  of  systems  of  this  type  with  periodic  coefficients  depending  on  a 
small  (dimensionless)  parameter,  e.g.. 


Al  (x/u),  Al  (x/u),  {xju) 


(for  0  <  1/  <<  l)  was  carried  out  in  (Bensoussan,  Lions,  and  Papanicolaou 
1978)  and  in  (Sanchez-Palencia  1980)  using  multi-scale  methods  (homog¬ 
enization).  The  periodic  case  represents  media  with  regularly  imbedded 


components. 


The  “inhomogeneities”  in  the  media  which  govern  the  scattering  of  EM 


radiation  from  foliage  covered  terrain  are  randomly  distributed,  and  so,  the 


homogenization  (multi-scale)  method  must  be  extended  to  treat  this  case. 
A  general  theory  for  partial  differential  equations  v/ith  random  coefficients 
was  developed  in  (Papanicolaou  and  Varadhan  1979).  This  work  forms  the 
basis  for  our  treatment  of  the  random  scattering  (and  absorption)  problem. 


3.2  Scaling  the  Model 


To  set  up  the  analysis,  we  must  prescribe  a  scaling  for  the  parameters 
of  the  random  medium,  in  this  case  the  permittivity,  permeability,  and 
conductivity.  The  general  form  is 

e{x)  =  e^{x)  =  e{x,x/i/) 

^i{x)  =  n''{x)  =  n{x,xlu)  (51) 

<7(1)  =  o‘'{x)  =  o[x,xlu) 

In  this  representation,  1/  >  0  is  a  small  (<<  1)  dimensionless  parameter 
which  chatracterizes  the  relative  scale  on  which  microscopic  variations  in  the 
medium  take  place.  In  previous  reports  we  have  shown  how  this  parame¬ 
ter  may  be  identified  in  terms  of  the  mean  free  path  between  elementary 
scatterers  in  a  foliage  covered  area.* 

The  presence  of  two  scales  in  (51)  suggests  that  we  aire  permitting  vari¬ 
ations  of  the  parameters  across  the  structure  of  two  different  spatial  types. 
The  “fast  scale”  y  =  x/v  captures  the  effects  of  microscopic  variations; 
that  is,  as  we  pass  from  one  elementary  scatterer  to  the  next.  Variations 
on  the  “slow  scale”  z  capture  the  effects  of  macroscopic  variations;  that 
is,  variations  as  we  move  from  one  type  of  vegetation  to  amother,  or  from 
one  region  of  a  given  type  to  another.  Note  that  as  y  changes  by  one  unit, 
*Th4t  if,  V  mmy  be  defined  in  terms  of  the  mean  free  path  between  scattering  cen¬ 
ters  normalised  by  a  characteristic  length  of  the  interaction  process,  e.g.,  the  length  of  a 
scattering  region  or  the  principal  axis  of  the  first  Fresnel  sone,  etc. 
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the  macroscopic  scale  z  changes  by  t/  <<  1  units;  and  as  z  changes  by  one 
unit,  y  changes  by  l/i/  >>  1  units. 

Using  this  scaling,  there  are  two  basic  cases  which  can  be  treated  by 
versions  of  the  method  of  multiple  scales: 

1.  The  functions  «(z,y),/x(z,y),<r(z,y)  are  periodic  in  their  second  ar¬ 
gument. 

2.  The  functions  €(z,y),;x(z,y),<y(z,y)  are  random  functions  of  their 
second  argument  which  are  stationary  and  satisfy  a  mixing  condition. 

The  first  case  was  treated  in  (Bensoussan,  Lions,  and  Papanicolaou 
1978)  and  in  (Sanchez-Palencia  1980),  (see  also  the  work  of  Tartar  and 
Murat  (1986))  among  many  other  related  problems.  The  second  case  pro¬ 
vides  a  basis  for  modeling  certain  types  of  foliage  covered  terrain. 

The  mixing  which  must  be  imposed  in  the  second  case  is  that  the  sta¬ 
tistical  correlations  between  random  events  in  the  medium  must  decrease 
rapidly  over  the  volume  of  interest.  That  is,  the  correlation  between  values 
of  «(z,yi)  and  €(z,y2)  decreases  as  jyj  —  y2l  — ►  oo.  In  effect,  this  means 
that  the  microstructure  of  the  foliage  at  one  point  yi  cannot  depend  too 
strongly  on  that  at  the  (distant)  point  yj.  This  is  clearly  reasonable  in 
many  situations.  It  may  not,  however,  be  the  case  in  wind  blown  grass. 

Notice  that  the  scattering  medium  is  still  permitted  to  have  (statistical 
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or  deterzninistic)  variations  which  take  place  in  the  “macroscopic”  x  spatial 
scale. 


3.3  The  Multiple  Scale  Hypothesis 


The  objective  of  the  multiple  scale  asymptotic  analysis  is  to  compute  ap¬ 
proximations  to  the  fields  E''{tyX)^H‘'{t^x)  in  the  limit  as  i/  — »  0,  that  is, 
in  the  limit  as  the  microscopic  variations  become  increasingly  dense. 

A  formal  procedure  for  computing  approximations  is  to  assume  an  ex¬ 
pansion  of  the  fields  in  the  form 

E''{x,xlu)  -  Eoix)  +uEi{x,y)  +  u^E2{x,y)  +  •••  (52) 

with  y  —  xfv  and  a  similar  expansion  for  H'' .  Introducing  the  change  of 


coordinates 


d  1  d 


dx  dx  V  dy 

in  Maxwell’s  equations,  substituting  the  expansions  for  E*'  and  iif*',  one 
finds  a  sequence  of  equations  for  the  terms  Ek{x,y),k  =  0,1,2,...  (and 
Hk)-  These  may  be  solved  using  a  device  called  a  “corrector”  which  is 
derived  from  a  kind  of  “separation  of  variables”  argument. 

If  the  dissipation  is  zero  [a  =  0),  the  result  is  a  representation  of  the 
solution  for  the  approximate  field  in  the  form 

d  f  9«(®)  0  1  [  ^o(x)  1  [  0  Vx  1  [  Eo{t,x)  1 


0  g^(i)  I  I  fro(i)  I  1  -Vx 


Ho[t,x) 


£;o(0,i)  =  £o(a:),-ffo(0,x)  =  Ho{x) 

V  .  [q,Eo]  =  0,  y  •  [g^Ho]  =  0 


The  (matrix-valued)  functions  9<(x),  ^^(x),  which  depend  on  the  macro¬ 
scopic  spatial  scale,  are  defined  in  terms  of  £(i, y), m(i, y)  by  the  “corrected 


averages 


q'J  (x)  =  M  a,y  6ij  -h 


dxfiy) 


5  J  —  2, 3 


where  onj  is  either  e.y  or  yLij,  and  the  function  x^iv)  is  the  corrector  associ¬ 
ated  with  £  or  fi.  Note  that  the  effective  parameters  in  the  approximation 
qt,q^  are  not  just  the  averages  of  the  rapid  variations  of  the  parameters 
£(x,*),At(x,')  over  the  medium.  They  include  the  “correctors”  which  re¬ 
tain  the  microscopic  interaction  effects  (multiple  scattering)  in  the  final 
approximation. 

The  correctors  satisfy  a  system  of  the  general  form  (see  the  following 
sections  for  details) 

V  •  [Q(ey  -t-  Vx®(y)]  =0,  M(xil  =  0  (56) 

where  ey  is  the  natural  basis  vector  in  In  these  equations  M  is  the 
operation  of  averaging  over  a  typical  cell  in  the  domain  for  the  periodic 
case;  and  it  is  expectation  with  respect  to  the  stationary  distribution  de¬ 
scribing  the  medium  in  the  random  case.  When  the  conductivity  <7^0, 
the  homogenized  model  is  more  complex,  including  a  “memory  effect"  -  see 
below. 


3.4  Implementation  of  the  Approximation 


Thus,  to  implement  the  approximation,  we  must  solve  the  equations  (56)  for 
the  correctors  and  then  compute  the  averages  in  (55).  Computation  of  the 
correctors  is  possible  in  the  periodic  case  (see,  e.g.,  (Begis,Duvaut,  Hassim 
198l)(Bougat  Derieux  1978)).  However,  in  the  random  case  of  interest  in 
this  project,  the  computation  is  very  difficult.  The  equations  (56)  in  the 
random  case  are  defined  path  by  path;  hence,  one  would  have  to  have  a 
complete  statistical  characterization  of  the  medium  to  be  able  to  solve  them. 
This  is  impossible  in  practice,  and  one  must  evaluate  the  correctors  by  a 
procedure  like  Monte  Carlo  simulations.  Hence,  these  equations  should  be 
regarded  as  the  basis  for  derivation  of  further  approximations  or  alternative 
representations  for  the  medium  which  require  only  physically  attainable 
statistical  information,  e.g.,  second  order  statistics  in  the  microstructure 
of  the  medium.  Such  approximations  have  been  derived  in  other  problems 
(conductivity,  porosity,  etc.)  in  the  past.  We  shall  consider  these  methods 
and  their  adaptation  to  the  EM  scattering/absorption  problem  in  the  next 
section. 

Before  giving  the  details  of  the  derivation  of  the  approximation,  we  shall 
make  a  few  remarks  on  its  interpretation. 

•  First,  note  that  the  general  form  of  the  approximation  is  valid  (with 
a  very  different  mathematical  interpretation)  in  both  the  random  and 


•  Second,  the  role  of  the  correctors  in  providing  the  “correct”  form  of 
the  approximation  cannot  be  dismissed.  Approximations  which  sim¬ 
ply  average  the  variations  of  the  permittivity  and  permeability  over  a 
domain  are  incorrect  -  making  an  error  of  “order  one”  in  the  small  pa¬ 
rameter  V.  That  is,  the  approximation  (55)  (56)  is  an  approximation 
of  “order  v”  (at  least  in  the  periodic  case);  omitting  the  correctors 
invalidates  this  estimate. 

3.5  Derivation  of  the  Effective  Parameters  Using  Ho¬ 
mogenization 


There  are  two  basic  cases  to  consider: 


1.  The  periodic  case  when  the  inhomogenieties  in  the  medium  are  spa¬ 
tially  periodic;  and 


2.  The  case  when  the  inhomogenieties  in  the  medium  are  r^mdomly  dis¬ 
tributed. 


The  expressions  for  the  effective  parameters  are  easier  to  grasp  in  the  peri¬ 
odic  case;  and  they  are  more  readily  computable.  For  this  reason  we  shall 
treat  this  case  first  as  background  to  the  random  case. 
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3.5.1  Periodic  Media 


Our  treatment  of  this  problem  is  based  on  (Sanchez-Palencia  1980).’ 

To  begin,  we  rewrite  Maxwell’s  equations 

dD'’ 

-^  =  VxH*'-J‘'  +  F  (57) 

at 

dB^ 

F*'(i,o)=o,  fr‘'(i,o)=o 

where  F,  G  are  (localized)  source  terms  and’ 

Dr  =  ««(j)£?.  =  j;  =,„(£)£;.  (ss) 

We  assume  that  c,  a  are  smooth  Y  periodic  functions  of  xju  which  are 
symmetric  and  positive  definite.^  Adopting  the  multiscale  hypothesis 

F‘'(i,t)  =  F°(i,t)  +  i/F‘(i,y,t)  H -  (59) 

f)  =  E\x,  t)  +  (i,  y,  0  +  •  •  • 

and  similarly  for  D*'  ,B''  and  J*'  with  each  term  in  the  asymptotic  expansion 
periodic  in  y.  We  shadl  use  the  change  of  coordinates 

dx  dx  V  dy 

^Se«  also  (Bsnsoussan,  Lions,  Papanicolaou  1978)  and  (Tartar  1979). 

*We  use  the  summation  convention  throughout  this  section  •  repeated  indices  in  an 

expression  are  summed. 

♦That  is,  «<y  =  tn,  >  'r|C|’,'T  >  0. 


with  y  =  x/i/,  and  similarly 


div  — »  div,  +  —  div„. 

u 


Rewriting  Maxwell’s  equations  in  terms  of  these  coordinates,  inserting  the 
expansions  (59)  into  them,  and  equating  coefficients  of  like  powers  of  u,  we 
arrive  at  a  sequence  of  equations  for  —  The  and 

0{t/°)  terms  are 

V,  x£r“  =  0,V,  xB“  =  0  (61) 

(62) 


^  =  +  VyX  -  r  +  F 


^  =  -V,xE°-V^xE'+G 

dt  *' 


The  equations  (61)  imply  that  Ho{x.,y,t),  Eo{x,y,t)  are  gradients  in  y  for 
fixed  values  of  the  other  arguments. 


The  terms  x  E^  and  x  are  derivatives  of  periodic  (in  y)  func¬ 
tions;  hence,  they  have  zero  mean  with  respect  to  the  averaging  operation 


■M(V,  X  ^‘)  =  jlq  jy  X  E'dy  =  0 


(63) 


(and  similarly  for  Vy  x  /f‘)  where  Y  is  a  typical  “cell”  of  the  periodic  struc¬ 
ture  and  \Y\  is  the  volume  of  the  cell.  Therefore,  averaging  the  system  (62) 
over  the  cell  Y,  we  arrive  at  the  homogenized  system 

dD^ 


dt 


=  V,x  +  F 


(64) 
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=  -V,  X  £?“  +  G 


where  the  overbar  denotes  the  averaging  operation.  To  complete  the  model 
we  must  derive  the  homogenized  constituent  equations  corresponding  to  (58). 
These  laws  take  two  different  forms,  depending  on  the  presence  or  absence 
of  dissipation  (a  =  0  or  not). 


Consider  the  divergence  of  (57) 


fdD''  \ 

div  =div  F 

{^)  =  ° 


Using  the  aruttaz  (59)  and  the  coordinates  (60),  we  have,  at  order 


which,  with  the  zero  initial  conditions,  implies 


div^B®  =  div„(/i(x,y)B°(i,y)]  =  0 


From  (61)  we  have  that  and  are  gradients;  hence,  we  csm  write 


E^-  E°  =  V^,  =  VV' 


with  the  means  of  ^  and  x(f  zero.  Using  these  in  (66)  (67),  we  have 


(71) 


These  are  the  key  equations  in  application  of  the  homogenization  method 
to  Maxwell’s  equations.  If  we  regard  (70)  as  an  equation  for  the  potential  0, 
then,  noting  that  does  not  depend  on  y,  we  can  solve  (70)  by  assuming 
the  “separation  of  variables” 

=  Xk{x,y,t)Ht{x,t)  (72) 

and  regarding  (70)  as  an  equation  for  x(y)  a  function  of  y,  treating  x,t  as 
fixed  parameters.  We  call  X*(^»y»0>^  =  lf2, ...  ,d  (d  is  the  dimension  of 
the  apace)  the  correctors  associated  with  Using  (70)  (72)  the  equation 
for  the  corrector  is 


This  problem  plays  a  key  role  in  the  homogenization  method.  We  call  it 
the  cell  problem.  The  assumption  that  /Xj,(i,j/)  is  periodic  in  y,  symmetric, 
and  positive  definite  (as  a  matrix  for  each  x,y),  guarantees  that  (73)  has  a 
unique  solution^  which  has  zero  mean 

The  solvability  condition  for  (73)  is  that  the  right  hand  side  have  zero  mean. 

This  holds  since  /i,y  is  periodic  and  smooth. 

*  Among  the  set  of  functions  periodic  in  y  with  Srst  derivatives  square  integrable  over  a 
cell,  i.e.,  periodic  functions  in  See  (Sanches-Palencia  1980),  pp.  61-54,  or  (Ben- 

soussan.  Lions,  Papanicolaou  1978). 
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The  homogenized  form  of  the  field  equations  may  be  obtained  using  the 
correctors.  Writing  out  the  0{u°)  term  in  the  expansion  of  (65),  we  have 


diVjB^  +  div,B“  =  div^G  (74) 

« 

Using  M  to  denote  the  averaging  operation,  the  solvability  condition  for 
this  equation  is 

M{div.B°  -  div,G}  =  0 


which  we  rewrite  as 

M  [div,{^(i,y)(Vy0  +  S°{x)\  -  G}]  =  0 
Using  the  representation  (72)  for  the  potential,  (75)  holds,  if 

>i[m(V,x*'  +  1)Ho-G]  =0. 

This  provides  the  definition  of  the  homogenized  magnetic  permeability, 
as 


(M.;(x,y) 


(76) 


The  effective  permittivity  and  conductivity  will  depend  on  the  frequency 
of  the  incident  radiation;  and  the  analysis  must  reflect  this.  Recall  (66), 
rewritten  as 


— €^o  +  oEq 


d 


)=0 
fv„<^  +  Eo 


i| 


(77) 


=  g^div,  +  div, 


(78) 
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r  a  .  ■)  r  _  1 

+div^  ^0  +  div^  [oEq^ 


Thus,  the  local  “cell”  equation  is  an  evolution  (in  t)  equation. 


To  treat  (78)  we  use  a  technique  introduced  in  (Sanchez-Palencia  Sanchez- 
Hubbert  1978).  Let  Vy  be  the  space  of  suitably  smooth  functions  fl(y)  which 
are  periodic  with  zero  mean,  and  let  Vy  be  endowed  with  the  inner  product 


f  d<f>  dO 
<  4>,6  >Vy=  / 

Jy  dyi  dyj 


Using  this,  (78)  becomes 


d  f  \  d<f>  /  &  \  1  . 

f  d<t>  dd  ,  f  /  s,  X 


Using  the  inner  product  notation,  this  may  be  written  compactly  as 


—  <  <^  +  cj(£'o)ji^  >Vy  +  <  >Vy  -t-  <  c^(£'o)y,^  >Vy  =  0  (80) 


where 


1  f  de 


z  f  de 

<  c;,fl  >Vy=  /  c 

Jy  dyi 


The  evolution  equation  (81)  can  be  solved  for  the  potential  <t> 


m  =  -c;(£„),(i)  +  jr‘»-''i‘-’ic;(£.)i(r)* 


cj  =  .4c;  -  c; 


Using  this  expression  for  we  can  compute  the  average  flux  vector  Do 
and  average  induced  current  Jo  a^d  define  the  effective  permittivity  and* 


conductivity.  Using  Eo  =  Eo  +  we  write  the  average  of  Do 

=  eEo  as 

(Do),  =  M  {£.,  (Fo  +  V^)} 

(82) 

+  +  /J  (r)rir) 

This  may  be  written,  compactly  as 

(Do),  =  <,(£„),(()  +  £  6:,.((  -  r)(£o)o(l  -  T)dr 

(83) 

where 

(84) 

In  a  similar  fashion,  we  can  compute  the  induced  current 

(85) 

(Jo),  =  <,(Fo)y(0  +  /J  b^,{t  -  r){Eo):{t  -  r)dr 

(86) 

where  a"  and  6*  are  defined  by  {85)(86)  with  e  replaced  by  a. 

Using  Laplace  transforms,  we  can  more  readily  express  the  effective 
scattering  parameters  as  a  function  of  frequency.  Indeed,  if  we  let  D-icj) 
be  the  Laplace  transform  of  D*'  with  respect  to  <,  and  similarly  for  the 
other  variables,  then  Maxwell’s  equations  take  the  form 

wl?*'  -I-  J"  =  V  X  H*'  -I-  F  (87) 
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The  claim  is  that  Maxwell’s  equations  approach  the  homogenized  laws 


ujD^  +  J^  =  (w£  +  (88) 

where  the  superscript  (•)*  denotes  the  “homogenized”  laws.  Using  the 
evolution  equations  (83)  (86)  we  can  identify 

uidij  +  cj6Jy(cj)  +  alj  +  Sfy(u;)  =  (we,;  +  >  0  (89) 

This  completes  the  derivation  of  the  (first  order)  homogenization  theory 
for  Maxwell’s  equations. 

Remark:  Using  the  arguments  in  (Bensoussan,  Lions,  Papanicolaou  1978) 
and  (Sanchez-PalenciaSanchez-Hubbert  1978)  and  (Sanchez-Palencia  1980), 
it  can  be  shown  that  the  homogenized  Maxwell’s  equations  (83)  (86)  (or  the 
frequency  domain  version)  have  a  unique  solution.  Moreover,  it  can  be 
shown  that  the  homogenized  Maxwell’s  equations  are  the  limits  as  i/  — »  0 
of  the  original  scaled  Maxwell’s  equations  (in  an  appropriate  weak*  topol¬ 
ogy). 

Rather  than  give  these  arguments,  we  shall  turn  our  attention  to  the 
case  when  the  coefficients  in  Maxwell’s  equations  are  random  processes  in 
the  spatial  vairiables. 
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3.5.2  Random  Media 


As  we  have  noted  earlier,  the  formulas  for  the  effective  parameters  in  the 
homogenized  form  of  Maxwell’s  equations  carry  over  in  form,  at  least,  to 
the  random  case.  The  averaging  operation  >1  is  interpreted  as  expectation 
with  respect  to  an  appropriate  probability  measure.  While  the  procedure 
of  simply  expanding  the  field  quantities  in  asymptotic  series  in  the  small 
parameter  is  only  formal  in  the  ramdom  case,  it  can  still  be  used  with  appro¬ 
priate  cautions  to  determine  equations  for  the  homogenized  field  quantities. 
The  arguments  for  convergence  of  the  scaled  field  quantities  to  the  homog¬ 
enized  quantities  are,  however,  substantially  different  in  the  random  case. 

The  central  difficulty  in  treating  random  coefficients  arises  in  the  mean¬ 
ing  and  evaluation  of  the  “correctors”  They  axe  random  pro¬ 

cesses  which  must  be  evaluated  (numerically)  by  a  Monte  Carlo  type  pro¬ 
cedure.  We  shall  take  up  this  issue  later  in  the  project.  In  this  section 
we  shall  examine  the  more  fundamental  issue  of  existence  and  uniqueness 
of  solutions  to  the  corrector  equations.  Since  each  of  the  three  correctors 
X'.X^.X"  satisfies  the  same  type  of  elliptic  equation  with  random  coeffi¬ 
cients,  we  shall  treat  the  generic  problem:  Find  x{y\'*>)  such  that 


x(0;u)  =  O.ExM  =0  Vy 


(80) 


^xl^U}  is  &  squaire  integrable  stationary  process  Vj  sind 


in  the  sense  of  distributions.  We  assume  that 

<  — lei'  ve6  3i",ao>o  (91) 

i,, 

where  a,j  (y,a;)  and  gj(y,uj]  are  square  integrable  stationary  processes,  t,  j  = 
1, . .  .n. 


Notation  -  Assumptions  We  shall  consider  a  set  up  for  stationary  pro¬ 
cesses  as  presented  in  Papanicolaou  and  Varadhan  (Papanicolaou  Vairadhan 
1979).  Let  (n,A,P)  be  a  probability  space  and  define  M  =  L^(n,Jt,P).  We 


assume  that 


M  is  separable 


there  exists  a  strongly  continuous  unitary  group  r„  on  K ,  y  €  9?"  (93) 

Ty  is  ergodic,  which  means  if  /  €  ^  satisfies  (94) 

Tp/  =  /,  Vy,  then  /  is  a  constant. 

>  0.  then  T^f  >  0  and  r„l  =  1  (95) 


The  group  Ty  has  a  spectral  resolution  defined  by 


r,  =  / 

Jit” 


i: 
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where  U{dX)  ia  a  projection  valued  measure.  We  consider  the  complex 
extension  of  )/,  provided  with  the  scalar  product  Efg,  whenever  /,  g  are 
two  elements  of  M.  The  measure  (7(A),  A  Borel  subset  of  SR”,  satisfies 

EU{^)fU{A')g  =0  v7,g€i7,AnA'  =  0 


and  by  ergodicity 


We  next  define 


EUiA)fU{A)g  =  EU{A)fg 


um)f  =  Ef 


a7m  =  ^  (r,7) 


which  axe  closed  densely  defined  linear  operators  with  domains  D  [Di)  in  )i 


Note  that 


E{gDi'f)  =  -E{D-gn  'i/rg  e  P(A), 


dJ  =  i  [  X^U{d\)f. 


If  Djf  =  0  Vy,  then  since 


=  f  X^E\U{dX)ff 

J9’' 

it  follows  that  U{dX)f  =0  VA  ft  0,  hence 

7  =  £^({o})7  =  E'f. 


Let  =  ny=i  ^i^j)  which  is  dense  in  M.  We  equip  with  the  Hilbert 
scalar  product 

=  ^/ff  +  E^DjD,g  (100) 

i=i 

We  identify  k  with  its  dual  and  call  k~^  the  dual  of  k^.  We  have  the 


inclusions 


k^  Ck  ck- 


each  space  being  dense  in  the  next  one  with  continuous  injection. 


The  family  is  also  a  strongly  continuous  unitary  group  on  k^,  since 

T,Dif  =  DiTj. 

It  can  be  extended  as  a  strongly  continuous  unitary  group  on  k~^  by  the 


formula 


<  Tj„ 7  >=<  7.,r_,7  >  V/  e  k\l  e  k- 


and  <,>  refers  to  the  duality  between  k^  and  k~^. 

Remark  1.  The  periodic  ease. 

Let  n  be  the  unit  n  dimensional  torus,  A  the  o-algebra  of  Lebesgue 
measurable  sets  and  P  Lebesgue  measure  on  fl.  Then  k  is  the  space  of 
measurable  periodic  functions  (period  1  in  each  component)  such  that 


/„(/(-))■ 


We  define 


du  <  oo. 


=  'f{'^  +  y) 


An  important  fact  in  the  periodic  case  which  does  not  carry  over  to  the 
random  case,  is  that  there  is  no  analogue  of  Poincare’s  inequality.  The 
consequence  is  following.  Consider  the  quotient  space  of  elements  of 

which  differ  by  a  constant.  Denote  by  [/)  the  equivalence  class  related 
to  an  element  /,  then  the  quotient  norm  is  given  by 


|7l  \\={E\'j-Eff  +  '£E\D,;]‘y>'‘ 


1/2 

This  is  not  equivalent  to  (2,  E\Djf\^'j  .  In  the  periodic  case  one  has 


1/1  1P=  /„  1/  - 


where  m(7)  =  f^  f{u]duj.  Poincare’s  inequality  implies  that  ||  [7]  ||  is 

1/2 

equivalent  to  (E,  In  ■ 


Consider  now  random  variables  not  necessarily  in  M .  We  assume  that 


Ty  is  a  linear  group  on  the  set  of  complex  random  variables  such  that 


•  •  -  Vk  complex  random  variables,  0  Borel  bounded  function  on  C*,  then 


E<t>{Tyfii,...,Ty^k)  =  E<f>{f}u...,fjt) 


(101) 


y,u>  Tyff  is  measurable, 


Tyi]  >  0  if  ^  >  0. 


A  stationary  process  is  a  stochastic  process  T](y;w)  which  csin  be  repre¬ 
sented  in  the  form 

T}{y\u)  =  Tyri{u).  (102) 

The  space  of  square  integrable  stationary  processes  can  be  identified 
with  U.  Moreover 

p-{y\u)  =  D,Tyri{u)  =  TyDif,  (103) 

oVi 

iffie 

Note  that  the  continuity  assumption  on  Ty  implies  that  the  square  inte¬ 
grable  stationary  processes  are  necessarily  continuous  function  of  y  with 
values  in  U .  Hence,  if  ^  €  )(,  ri{y',ai)  €  C°(IR";  )/),  space  of  uniformly  contin¬ 
uous  functions  on  SR"  with  values  in  ){.  If  fj  6  then  Tj{y;ui)  6  C'(SR";  U). 

Note  that  Di  G  If  r)  G  we  can  consider  the  distribution 

derivative  with  values  in  )( ,  defined  by 

/..  = - L  ^ 

Let  us  check  that 

-  QQ  m 

-  rj{y,u)—dy  =  j  DiTyfi9{y)dy  (104) 

y*"  ayi  Jm- 

which  is  an  equality  in  This  proves  that  r7(y,i4;)  G  C“'(SR";  )/“*)  and 
the  distribution  derivative  dyjdyi  with  values  in  U  can  be  considered  as  a 
continuous  function  with  values  in 
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■v  .“S  -V 


To  prove  (104)  pick  ^  then 

[  D,Ty,^e{y)dy  =  -E  f  Di^Tyf]e{y)dy 
Jit*  Jit* 

=  -E  f  DiT-yfiB{y)dy 


=  '^{~y,^)vO{y)dy 

f  90 

=  -E  <<>(-y,u;)T7  — {y)<iy 

oyi 

f  -  9$ 

=  -E  J^^^T]{y,u)  —  {y)dy 


which  completes  the  proof  of  (104). 


The  Cell  Problem  We  consider  here  stationary  processes  a,;  (y;w)  such 
that 


ool^l’ <  <  -i^l*  Ve6lR",ao>0  (105) 

Let  y;(y,  w)  =  Tyj,  square  integrable  stationary  processes,  j  =  1, . . .  n.  We 
shall  solve  the  problem:  Find 


X(y;w)  €  C^(!R'‘;»),x(0;u;)  =  0,£x(y)  =0  Vy 
^xl^Vi  i«  *  square  integrable  stationary  process  Vj 


(106) 
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in  the  sense  of  distributions  with  values  in  ){  (or  as  continuous  functions 
with  values  in 


Papanicolaou  and  Varadhan  (Papanicolaou  Varadhan  1979)  have  shown 
the  existence  and  uniqueness  of  the  solution  of  (106).  We  shall  reproduce 
their  proof  with  minor  changes.  Note  that  x(y‘»‘*')  itself  is  not  a  stationary 
process.  This  is  a  big  difference  with  respect  to  the  periodic  case  and  relates 
to  Remark  1.  Note  also  that  xCVit*')  ^  C7*(5R";  ^~*).® 

Let  Xj  €  H  such  that 

we  can  assert  that 

=  ExkDj^,  6  )i^.  (107) 

Indeed  we  have  to  show  that 

DkXi  =  EfjXk  (108) 

as  an  equality  in  But 

TyDkXj  = 

5*X  /  ^ 

=  o  a~  (Vi 
dyjOVk 

hence 

'EyDkXj  ~  "EyDjXk 

*By  virtue  of  (95)  it  ia  sufficient  to  hsve  (106)  for  y  =  0. 
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which  implies  (108). 


hence,  as  is  etisily  seen 


Because  of  (108)  we  have 


EoijXjDi^^  =  0. 


(Ill) 


f  XkU{dX)xi  ^  f  XjU[dX)xk 


which  implies 


-04>^=  I  E 


(-tA,-^)(tAt-^) 

1*-A  -  0\^ 


U{dX)x,{u^) 


(-tA,- -/?)(, A,-/?) 


ItA  -  0\^ 


U{dX)xk{<^)  =  Xk- 


Therefore,  (111)  reads 


EdyXiXi  +  PEayX;^^  =  0. 


(112) 


However, 


j.k  I*''  -  ^1 

-+  E^^(0}X;Xt  as  0  ^0 

},k 

and  by  ergodicity  and  property  (109),  we  get 


Therefore, (112)  implies 


-  Oa30  -  0. 


EdijXiXi  =  0 


hence  x,-  =  0,  which  implies  also  x  =  0. 


Existence: 


Let  /3  >  0,  we  solve  the  problem 


^  5yo  dyj 

X^{y;u;)  =  TyX^,x^  e 


This  problem  is  equivalent  to 

EoiyDjX^ Di^  +  0Ex^^  =  —EgjDj^y  6 


We  easily  deduce  the  estimates 

E\D^X^\^  <  C 
PE{xy  <  C. 

Let  \i8  extract  a  subsequence  such  that 

DjX^  — »  X}  H  weakly. 

EDjX^Dki)  =  EDkX^Dj^ 

EDjx''  =  0 

we  deduce  (108)  and  (109).  Going  to  the  limit  in  (114)  we  have 


Define  then 


xM  =  [  {«■'"- (lie) 

then 

=  ^»Xk 

and  x(0;‘*/)  =  Oi  ■E'x(y)  =  0-  Then  (115)  can  be  written  as 

-D,a,yXj  =  DiQi  equality  in 
which  is  indeed  (106). 


4  Bounds  for  Effective  Parameters 

In  previous  sections  we  have  examined  methods  for  representation  and  com¬ 
putation  of  effective  parameter  models  for  heterogeneous  media  in  terms 
of  statistical  characterizations  of  the  media.  These  methods  lead  to  repre¬ 
sentations  for  the  effective  parameters,  e.g.,  the  effective  dielectric  constant 
and  conductivity,  which  are  difficult  to  evaluation  in  any  but  the  simplest 
cases.  In  this  section  we  examine  the  derivation  of  bounds  on  the  effective 
parameters  as  a  means  of  quickly  providing  good  approximations  to  the 
behavior  of  the  medium. 

There  has  been  a  great  deal  of  work  in  this  area,  not  only  in  electromag¬ 
netics,  but  also  in  conductivity  and  elasticity.  This  work  dates  from  the 
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early  years  of  this  century  (Bergman  1978),  and  it  includes  some  very  recent 
work  (e.g.,  (Ericksen,  et  al.  1986)).  We  have  been  studing  the  portions  of 
this  work  applicable  to  the  case  at  hand  -  scattering  and  absorption  of  EM 
radiation  by  foliage  covered  terrain.  Several  key  issues  must  be  resolved 
before  the  work  can  be  applied  directly  to  scattering/absorption  from  fo¬ 
liage  with  more  than  one  kind  of  elementary  scatterer.  We  shall  point  out 
these  issues  in  the  following  paragraphs. 

4.1  Definition  of  the  Effective  Parameters 

The  starting  point  is  a  definition  of  the  effective  parameter,  e.g.,  the  effec¬ 
tive  permittivity.  Suppose 

€{x,uj]  =  €ili(u;)  -t-  •  •  •  -I-  e/v'lvit*')  (lif) 

is  the  dielectric  “constant”  of  a  composite  medium,  where  €j,i  =  1,. . .  ,.V 
is  the  (complex)  dielectric  constant  of  medium  t  ,  and  the  indicator  function 
1,  is  1  for  all  samples  w  €  H  which  have  medium  t  at  point  z,  and  zero 
otherwise.  We  use  (n,/,P)  to  denote  the  probability  space  on  which  the 
statistics  of  the  phenomena  are  defined.  The  quantity 

fj  =  j^P{du;)lj{u)  (118) 

is  the  volume  fraction  occupied  by  the  material  of  type  j  =  1,2,. ..  ,N. 

The  spatially  dependent  random  conductivity  and  magnetic  permeabil¬ 
ity  are  similarly  defined.  In  our  case  these  functions  depend  on  the  random 
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geometry  of  the  scattering  foliage  which  may  be  available  only  in  terms  of 
second  order  statistics.  The  bounds  obtained  for  the  “effective  parameter” 
approximations  of  these  functions  should  not  depend  on  statistics  beyond 
this  level. 

If  E{x^u))  and  D{x,<jj)  are  the  electric  and  displacement  field  vectors 
(stationary  in  time),  then 

A(x,w)  =  (119) 

V  •  D{x,uj)  =  0 

V  X  E{x,ijj)  =  0 

and 

I  P{<Lj)E{x,uj)  =  E 
Jn 

provides  the  boundary  condition,  where  .E  is  a  constant  field  incident  on  the 
region.  By  normalizing  the  field  magnitudes  and  redefining  the  coordinates, 
we  can  take  E  =  Ck,  the  unit  vector  in  SR^.  Let  E*  and  E*  denote  the 
corresponding  fields.  Then  (119)  may  be  rewritten 

Di{x,u;)  =  ^  £„(i,u;)E*(x,u;)  (120) 

0=1 

The  effective  dielectric  constant  may  be  defined  as 

=  (121) 

In  (Golden  Papanicolaou  1983)  it  was  shown  that  this  essemble  average 
coincides  with  the  usual  definition  involving  a  volume  average. 


Equation  (121)  may  also  be  written  in  the  symmetric  form 

6**  =  /  P{du)'tE^Ey  (122) 

where  the  bar  denotes  complex  conjugate.  This  expression  suggests  the 
interpretation  of  e*  as  the  dielectric  constant  of  a  fictitious,  homogeneous 
medium  (for  the  volume  V)  which  provides  the  same  value  for  the  electro¬ 
static  energy  stored  in  the  volume  occupied  by  the  heterogeneous  medium.^ 

4.2  Methods  for  Computing  Bounds 

Approximation  to  the  effective  permittivity  (or  conductivity  in  thermal 
systems)  have  been  derived  by  three  basic  methods.  The  original  work  of 
Haskin  and  Shtrikman  (1962)  was  based  on  variational  principles,  and  these 
have  been  reformulated  in  recent  years  to  provide  one  of  the  basic  meth¬ 
ods  for  deriving  aproximations  (Milton  1981a, 198 lb).  The  second  method 
for  deriving  bounds  is  the  method  of  “compensated  compactness”  (Tartar 
1986).  The  third  method  involves  the  use  of  complex  analysis  and  “rep¬ 
resentation  formulas”  (Bergman  1978)  (Golden  1986)  (Golden  Papanicolaou 
1983).  We  shall  review  the  first  two  methods  in  detail  in  another  report. 
Here  we  shall  briefly  describe  the  analytical  approach  to  deriving  bounds, 
indicate  its  range  of  applicability,  and  discuss  our  ideas  for  adapting  the 

’^The  volume  average  form  ia  =  [Jy  e(r)lE(r)|’(iV]  /mea»(V)  where  V  ia  the 

repreaentative  region,  mea»(V)  ia  ita  volume,  and  Eo  =  [fy  E{r)dV]  /meas{V)  ia  the 
average  (incident)  field. 


method  to  treat  scattering  and  absorption  from  foliage. 


If  we  substitute  (117)  into  (121)  and  divide  by  e^,  then  we  obtain 


. h^_i)  =  J  E^iuj)  (123) 


Here  h,  =  u/es-  The  key  idea  in  the  work  of  Bergman  (Bergman  1978)  is 
to  regard  m(-)  as  a  function  of  the  complex  variables  hi, . . . ,  —  1  in  C^~^. 

This  permits  one  to  use  (spectral)  representation  formulas  for  m(-)  which 
allow  the  derivation  of  bounds  based  on  the  extreme  points  of  certain  sets 
of  measures.  By  expanding  the  representation  formulas  about  the  case  of  a 
“homogeneous”  medium  hj  —  \,j  =  1, . . . ,  iV  —  1,  it  is  possible  to  obtain  a 
characterization  of  the  underlying  measures  in  terms  of  their  “moments.” 
These  moments  are  functions  of  increasing  amounts  of  information  on  the 
statistical  properties  of  the  medium.  For  example,  if  only  the  volume  frac¬ 
tions  /;  from  (118)  are  known,  then  only  the  first  term  in  the  representation 
can  be  computed  explicitly.  This  leads  to  the  classical  Wiener  bounds: 


For  two-component  media 


<  f*  £  /ifi  +  hf-i 


(124) 


For  three-component  media  - 


—  ^  <  c*  ^  /i«i  +  /jCj  -I-  /3C3 

«1  €J  «3/ 


(125) 


These  bounds  are  achieved  by  parallel  plane  configurations  of  the  materials. 


If  additional  information  is  available  on  the  statistical  description  of  the 
medium,  then  more  sophisticated  bounds  can  be  obtained.  For  example,  for 
a  two-component  medium  which  is  statistically  isotropic  the  first  moment 
of  the  spectral  measure  can  be  computed  (in  the  expansion).  This  leads  to 
the  Hashin-Shtrikman  bounds  (assume  without  loss  of  generality  Ci  <  fz) 

“  i/(<.  -  2 


4.3  Multicomponent  Media 


There  is  a  substantial  difficulty  in  passing  from  two-component  materials  to 
iV-component  materials  with  N  >  3  when  using  the  analytical  continuation 
method.  Basically,  one  must  deal  with  functions  of  more  than  one  complex 
variable  (hi,  hj, . . .).  The  precise  difficulty  arises  from  the  fact  that  the 
extreme  points  of  the  associated  set  of  spectral  measures  axe  not  known. 
Prior  to  the  paper  (Golden,  Papanicolaou  1985)  there  was  no  systematic 
method  for  treatment  of  multi-component  media  (using  the  continuation 
method).  In  the  case  of  three-component  media,  it  is  possible  to  circum¬ 
vent  this  problem  by  imposing  a  linear  relationship  among  the  complex 
parameters,  and  expanding  the  spectral  measure  about  the  (linear)  param¬ 
eter.  In  (Golden,  Papanicolaou  1985)  an  unconfirmed  technical  hypothesis 
was  introduced  which  permitted  recovery  of  the  Hashin-Shtrikman  bounds 
for  three-component  media.  Assuming  ci  <  €2  <  ts  we  have 
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These  bounds  are  “optimal”  for  materials  consisting  a  mixture  of  spheres  of 
all  different  sizes  of  ti  and  (2  materials  each  coated  with  £3  material  in  the 
appropriate  volume  fraction.  Additional  bounds  were  obtained  in  (Golden 
1986)  for  complex  permittivities  (thus  including  the  absorption  of  energy 
in  the  material). 

The  analytical  continuation  method  is  a  powerful  technique  for  deriva¬ 
tion  of  bounds  on  the  material  properties  of  multi-component  media;  and  it 
is  useful  to  examine  its  potential  for  the  evaluation  of  scattering/absorption 
phenomena  in  the  interaction  of  EM  radiation  and  foliage.  Our  evaluation 
of  this  methodology  is  not  yet  complete;  however,  the  following  points  are 
worth  making  at  this  time: 

•  The  method  is  based  on  the  assumption  that  the  medium  is  nearly 
homogeneous;  that  is,  the  normalized  permittivities  tj/e.v  ~  1-  This 
is  a  restrictive  assumption  for  foliage,  and  it  may  limit  application 
of  the  approximation  to  certain  kinds  of  foliage  (mixtures  of  grasses, 
etc.). 

•  While  the  case  of  two  component  media  is  reasonably  well  developed; 
the  case  of  three-component  media  is  less  clear.  This  is  the  key  case 
if  one  is  to  understand  the  interaction  of  scattering  effects  in  an  envi¬ 
ronment  of  two  different  types  of  scatterers  in  a  uniform  background. 
The  limitations  on  the  method  are  both  technical  and  physical. 


•  The  technical  limitations  are  mainly  due  to  the  incomplete  character¬ 
ization  of  the  set  of  spectral  measures  for  functions  of  two  complex 
variables. 

•  The  physical  limitations  are  more  fundamental.  The  method  does 
not  account  for  the  geometry  of  multiple  scattering  processes.  That 
is,  while  it  is  capable,  in  principle,  of  including  higher  order  sta¬ 
tistical  functions  in  the  expansion  for  the  elfective  permittivity,  the 
method  develops  the  bounds  in  a  fashion  which  is  well  removed  from 
the  physical  properties  of  the  medium.  The  extreme  points  of  the 
set  of  measures  which  are  used  to  determine  the  bounds  are  singular 
measures  (“delta  functions”)  which  may  not  correspond  well  to  the 
distributions  of  scatterers  in  typical  foliage.  It  would  be  better  to 
develop  the  bounds  in  a  (restricted)  set  of  meEisures  which  are  repre¬ 
sentative  of  foliage.  It  appears  to  be  possible  to  adapt  the  analytical 
continuation  method  to  this  case. 

We  are  continuing  our  investigation  of  this  methodology  for  the  con¬ 
struction  of  bounds  with  an  eye  to  improving  on  these  limitations.  Our 
objective  is  to  compare  the  approximations  produced  by  this  method  with 
those  derived  through  the  other  techniques  we  have  studied  and  develop 
a  hierarchy  of  approximations  for  scattering  in  foliage  (of  various  types). 
We  shall  design  a  series  of  numerical  experiments  to  test  the  methods  on 
representative  data. 
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